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ABSTRACT 

Context. Cyg OB2 is one of the most massive associations of O-type stars in our Galaxy. Despite the large interstellar reddening 
towards CygOB2, many studies, spanning a wide range of wavelengths, have been conducted to more clearly understand this as- 
sociation. X-ray observations provide a powerful tool to overcome the effect of interstellar absorption and study the most energetic 
processes associated with the stars in Cyg OB2. 

Aims. We analyse XMM-Newton data to investigate the X-ray and UV properties of massive O-type stars as well as low-mass pre- 
main sequence stars in Cyg OB2. 

Methods. We obtained six XMM-Newton observations of the core of Cyg OB2. In our analysis, we pay particular attention to the 
variability of the X-ray bright OB stars, especially the luminous blue variable candidate Cyg OB2 #12. 

Results. We find that X-ray variability is quite common among the stars in CygOB2. Whilst short-term variations are restricted 
mostly to low-mass pre-main sequence stars, one third of the OB stars display long-term variations. The X-ray flux of Cyg OB2 #12 
varies by 37%, over timescales from days to years, whilst its mean log j^- amounts to -6.10. 

Mjol 

Conclusions. These properties suggest that CygOB2 #12 is either an interacting-wind system or displays a magnetically confined 
wind. Two other X-ray bright O-type stars (MT91 516 and CPR2002 Al 1) display variations that suggest they are interacting wind 
binary systems. 

Key words. Open clusters and associations: individual: CygOB2 - Stars: early-type - stars: pre-main sequence - stars: individual: 
Cyg OB2 #12 - X-rays: stars 



1. Introduction 

The CygnusOB2 association is one of the richest and most 
massive OB associations known in our Galaxy. These prop- 
erties make it an interesting target for observations over a 
broad range of wavelengths from the radio to the y-ray do- 
main. Unfortunately, photometric and spectroscopic studies of 
Cyg OB2 in the UV (e.g. Herrero et al. 2001) and optical wave- 
band (e.g. Massey & Thompson 119911 Torres-Dodgen et al. 
119911 and Hanson 2003 ) are seriously hampered by the heavy 
and patchy absorption by dense molecular clouds. Nevertheless, 
the interest in this association has been boosted over the past 
decade by near-infrared studies, which have obtained a far more 
complete census of its stellar population, and by the discovery 
of two y-ray sources, a GeV and a TeV source, in this direc- 
tion. From 2MASS near-IR data, Knodlseder (120001 concluded 
that Cyg OB2 probably harbours about 2600 + 400 OB stars, 
among which 120 ± 20 are O-stars, within an angular radius 
of about one degree. Although more recent studies suggest a 
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slightly lower number of 90 - 100 O-type stars (Comeron et 
al. 12002). CygOB2 is certainly quite rich in massive stars and 
contains stars as early as 03 If* (Cyg OB2 #7). Accordingly, it 
has been suggested that the stellar winds of this concentration 
of massive stars in Cyg OB2 could be the origin of the slightly 
extended, unidentified TeV source TeV J2032+4130 (Aharonian 
et al. 120051 Albert et al. 120081) . although alternative scenarios 
have also been proposed (see Butt et al. [2008 and references 
therein). In the GeV energy domain, the unidentified EGRET y- 
ray source 3EG J2033+41 18 had also been proposed to be as- 
sociated with interacting winds from massive stars in the core 
of the association (Romero et al. 1999 Rauw 2004), but recent 
Fermi-L AT data have identified this object as a young energetic 
pulsar (PSR J2032+4127, Abdo et al. |2TJT0l . 

Four CygOB2 members have been extensively studied in 
the radio domain. Three of them are non-thermal radio emit- 
ters (Cyg OB2 #5, #8a, and #9, see e.g. De Becker 120071 and 
references therein) and are now known as binary or multiple 
systems hosting wind interaction regions (Kennedy et al. 2010 
De Becker et al. 120041 Naze et al. 120081 ). whilst the fourth one 
(Cyg OB2 #12) is a luminous blue variable (LBV) candidate and 
could be one of the most luminous and optically brightest stars 
in our Galaxy (e.g. Massey & Thompson 1 199 II) . 

X-ray emission from a handful of massive stars in Cyg OB2 
was discovered serendipitously when the EINSTEIN observa- 
tory was pointed at Cyg X-3 (Harnden et al. Q979 ). All sub- 
sequent major X-ray satellites have observed this association. 
The first CCD X-ray spectra of the four brightest massive stars 
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(CygOB2 #5, 8a, 9, and 12) in the association were obtained 
with ASCA-SIS (Kitamoto & Mukai [19961 The Chandra ob- 
servatory has observed CygOB2 several times, both with the 
HETG, to gather high-resolution X-ray spectroscopy of the 
brightest sources (Waldron et al. 2004), and the ACIS-I instru- 
ment, to collect CCD images and spectra of the members of the 
associ ation ( Albacete-Colombo et al. I2007al I2007bl Wright & 
Drake [20091 Wright et al. 2010 Q Thanks to the very narrow 
PSF and the low background, the lOOksec ACIS-I exposure of 
the central region of Cyg OB2 (centered near Cyg OB2 #9) re- 
vealed about 1300 sources (Wright & Drake 12009b . Most of 
these sources are very faint and are likely associated with pre- 
main sequence stars. This is unsurprising as IRAS and 12 CO 
surveys revealed the Cyg OB2 region to be a highly active star- 
forming region (Odenwald & Schwartz 119931 Parthasarathy et 
al. [19921 

The distance to the CygOB2 association remains some- 
what controversial. Massey & Thompson (1991) derived a dis- 
tance modulus of 1 1.2 ± 0.1, in agreement with the more recent 
mean spectro-photometric distance obtained by Kiminki et al. 
(2007). Hanson (120031) used a cooler O-star temperature scale 
than Massey & Thompson d 199 1 b and inferred a lower distance 
modulus of 10.4 mag (implying luminosities that are lower by 
a factor 2). An even lower distance modulus (9.8 mag) was ob- 
tained from the analysis of the light curve of the eclipsing bi- 
nary Cyg OB2 #5 (Linder et al. 12009b . although this latter result 
might be biased by the presence of third light as Cyg OB2 #5 
most likely hosts four stars (Kennedy et al. 120 101 ). Throughout 
this paper, we adopt the 10.4 mag distance-modulus of Hanson 
( 2003 ), corresponding to a distance of 1 .2 kpcj^jalthough we em- 
phasize that most of our conclusions are essentially independent 
of the distance. 

We report the results of our analysis of six XMM-Newton 
observations of the core of Cyg OB2. The main purpose of this 
campaign was to study the X-ray properties and variability of 
the X-ray bright massive members of the association. Large X- 
ray luminosities, X-ray variability and/or spectral hardness are 
frequently considered as indicators of either interacting wind bi- 
naries or single stars with magnetic fields (see e.g. Gagne et al. 
1201 11 ). although this is certainly not a one to one relationship 
and some O-star binaries actually display X-ray properties quite 
similar to those of normal single O-stars (e.g. Naze et al. 1201 11 ). 
When interpreting our data, we thus have to keep in mind the 
latest results of multiplicity studies in Cyg OB2 (Kiminki et al. 
12009 1 and references therein). 

This paper is outlined as follows. In Sect.[2j we present our 
data set and the data reduction. The general properties of the X- 
ray and UV sources in our field of view as well as their long and 
short term variability are discussed in Sect. [3] Section[4]presents 
a discussion of the X-ray emission of five relatively bright OB 
stars and Sect.[5]summarizes our conclusions. 

2. XMM-Newton observations of Cyg OB2 

Six pointings centred on CygOB2 #8a were taken with the 
XMM-Newton satellite (Jansen et al. 1200 II ). The first four obser- 
vations were obtained in October - November 2004 with roughly 
ten days between two consecutive pointings. Two follow-up ob- 

1 Most recently Cyg OB2 has been the target of a Chandra legacy 
survey (Drake 120 111 . 

2 This value agrees with the distance estimate of 1.32+^ kpc to- 
wards the W75N star forming region (Rygl et al. I2010i . which is part 
of the Cygnus X giant molecular cloud, as is Cyg OB 2. 



servations were gathered in April - May 2007 (see Table[T|). The 
EPIC cameras were used with the medium filter to reject optical 
light. 

The raw data were processed with SAS software version 
10.0. Several observations (especially the last two) were af- 
fected by high background events (so-called soft-proton flares). 
These episodes were rejected in our processing and the remain- 
ing total useful exposure times were 105ksec for both EPIC- 
MOS (Turner et al.[200TT i detectors and 74ksec for the EPIC-pn 
(Striider et al. 2001). Images were extracted for the soft (0.5 - 
l.OkeV), medium (1.0 - 2.0keV), and hard (2.0 - 8.0keV) en- 
ergy bands. A few stray-light features, caused by singly reflected 
photons from CygX-3, affect a small region in the lower right 
corner of the images. These parts of the images were excluded 
from the subsequent analysis (see Fig.[T|. 

We note that the best-fit positions of Te V J203 2+4130 are 
just outside our field of view (Aharonian et al. 2005, Albert et al. 
2008 ). Apparently diffuse X-ray emission from this TeV source 
was discovered with another XMM-Newton observation (Horns 
etal. 120071) . 

Some preliminary results for the first four XMM-Newton 
pointings were discussed by Rauw et al. (2005), whilst full anal- 
yses of the data for Cyg OB2 #8a, C ygOB2 #5, and Cyg OB2J9 
were presented by De Becker et al. (12006b . Linder et al. ([2009 ), 
Naze et al. ( 12010b , and Blomme et al. ( 12010b . These three objects 
are not discussed any further in this paper. 




Fig. 1. Energy-coded three-colour image built from our six 
XMM-Newton observations of CygOB2. The red, green, and 
blue colours correspond to the soft, medium, and hard energy 
bands used throughout this paper (see text). The individual EPIC 
images were exposure corrected, the stray-light features from 
Cyg X-3 (lower right corner) were excluded, and the out-of-time 
events from the bright central source (Cyg OB2 #8 a) were re- 
moved before the images were combined. A preliminary version 
of this figure was published by Rauw et al. ([2008 ). 

The Optical Monitor (OM, Mason et al. 12001b instrument 
onboard XMM-Newton observed Cyg OB2 through the UV Wl 
and UVM2 filters (see Fig. [2]). The first four pointings were per- 
formed in 'full frame low-resolution' mode, whilst the last two 
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Table 1. Journal of the six XMM-Newton observations discussed in this paper. 



Obs# 


Date 


Exposure time 


Clean time 


OM exposure time 




JD-2450000 


(ksec) 


(ksec) 




(ksec) 






MOS 


pn 


MOS 


pn 


UVW1 UVM2 


1 


3308.579 


19.7 


18.0 


19.1 


15.4 


3.2 


8.0 


2 


3318.558 


21.7 


20.0 


21.5 


17.6 


4.2 


9.0 


3 


3328.543 


23.7 


22.0 


23.0 


19.3 


5.0 


10.0 


4 


3338.506 


21.7 


20.0 


13.7 


10.9 


0.0 


8.2 


5 


4220.355 


30.5 


28.9 


11.8 


3.6 


10.3 


12.3 


6 


4224.170 


31.6 


30.0 


15.9 


7.0 


9.9 


13.7 



Notes. Columns 5 and 6 provide the useful exposure times after discarding time intervals of high background. 



Cyg 0B2 : OM data in UVW1 filter 



Cyg OB2 : OM data in UVM2 filter 
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Fig. 2. Mosaic image of the core of CygOB2 as seen with the Optical Monitor through the UVW1 (left) and UVM2 (right) filters 
during observation 5. The impact of the interstellar reddening on the number of sources detected in the UVM2 band is clearly seen. 



used the 'imaging mode' consisting of a mosaic of five expo- 
sures to cover the full 17 arcmin square field of view (FoV) of the 
instrument (Fig. [2]). The data were processed with the relevant 
commands of the SAS software. The OM images display some 
artifacts such as smoke-ring ghosts caused by light reflected off 
a chamfer in the detector window and modulo-8 pattern around 
bright sources (see Mason et al. 2001). The OM detector has 
a read-out time of 1 1 ms. The brightest sources in our field of 
view therefore suffer from coincidence losses that amount typi- 
cally to 10% for sources with 10 counts s (Mason et al. 1200 11 1. 
These losses are partially recovered during the data process- 
ing to provide 'corrected' count rates, which are then converted 
into AB magnitudes (Oke |1974| l. In our cases, a count rate of 
10 counts s _1 corresponds to UVWl and UVM2 AB-magnitudes 
of 16.066 and 14.912. The astrometry of a dozen O-type stars, 
spread all over the OM field of view, was compared to the coor- 
dinates taken from the catalogue of Massey & Thompson (1991 
hereafter MT91). We found a small shift (around 1.3arcsec in 
right ascension and 1.8 arcsec in declination) for each observa- 
tion. There was no obvious trend in this shift with position across 
the FoV, thus we applied a simple translation to the OM coordi- 
nates before performing cross-correlations with optical or X-ray 
catalogues. After this correction, the average distance between 
the OM sources and their counterpart in the MT91 catalogue 
amounts to 0.63 arcsec for 225 objects (see Sect. [3]l. 



3. The X-ray sources: general properties 

A three colour X-ray image derived from our combined data set 
is shown in Fig.[T] It can be seen that the faint X-ray sources 
are found all over the field of view, although some areas contain 
larger concentrations of sources. Bica et al. (2003) reported the 
discovery of two open clusters in the core of Cyg OB2 roughly 
centered on stars #8a and #22, respectively. Whilst we find a rel- 
atively large number of X-ray sources within the regions defind 
by Bica et al. (2003), we note that there are other places inside 
our field of view with similar concentrations of X-ray sources 
(such as around 20h32m47s +41°24'45"). 



3. 1. The nature of the X-ray sources 

The SAS routines detected a total of 206 sources in the combined 
EPIC images (outside the region affected by the straylight) with 
a significance threshold of 2Qjand average count rates as low as 
0.0004 and 0.0002 cts s 1 for EPIC-pn and EPIC-MOS. From a 
histogram of the number of detections versus count rate, we con- 
clude that our detections are complete down to rates of 0.0022 



3 This implies a probability of < e~ 20 that a random Poissonian fluctu- 
ation could have caused the observed source counts inside the detection 
cell. Over the entire EPIC field of view, the expected number of spu- 
rious detections amounts to about 2% of the total number of detected 
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and 0.0008 cts s _I for EPIC-pn and EPIC-MOS, respectively. All 
these sources were verified by eye and 199 of them were even- 
tually retained as real sources. There are a number of additional, 
fainter sources that were not selected by the detection algorithm, 
because they are located in crowded regions and their signifi- 
cance over the background therefore falls below our threshold. 
Finally, we stress that a few sources in our list could actually be 
blends of two fainter, spatially unresolved sources. Table[2] lists 
the coordinates of our sources, their X-ray properties, and the 
results of cross-correlation with various catalogues (see below). 
The full version of this table is available at the CDS. 

We cross-correlated our final list of X-ray sources with the 
2MASS point source catalogue (Cutri et al. 2003 Skrutskie 
et al. l2.0061 >. On the basis of the technique of Jeffries et al. 
( 11997b . the optimal cross-correlation radius, which offers the 
optimal compromise between obtaining the maximum number 
of true identifications and avoiding contamination by spurious 
coincidences, was established to be 4arcsec. For this correla- 
tion radius, we expect to find fewer than 5% of false cross- 
identifications. We also correlated our list of sources with the 
U BV photometric catalogue of Massey & Thompson dl99U 
and the list of OB star candidates provided by Comeron et 
al. (120021 . Twenty-four positive matches were found with the 
MT91 catalogue, to which we can add CygOB2 #5, which is 
not included in the compilation of Massey & Thompson dl991| l. 
Three OB star candidates listed by Comeron et al. (2002J fall 
within the field of view of our data. Two of them, CPR2002 A 1 10 
and A43 (= Cyg OB2 #16 = MT91 299), are likely counterparts 
to X-ray sources, whilst the third (CPR2002A17) is not de- 
tected. Furthermore, among the four Bry emitters (CPR2002 B5, 
B8, Bll, B13, in addition to Cy g OB 2 #9) in our field of view 
listed by Comeron et al. (2002), only star B8 is detected, al- 
though at a rather large angular distance from the position of the 
X-ray source. We find that a total of 20 X-ray sources correspond 
to O-type stars, and 5 are associated with B-type stars, includ- 
ing the B supergiant and LBV candidate Cyg OB2 #12. In addi- 
tion to Cyg OB2 #5, #8a, and #9, there are 5 OB stars (Cyg OB2 
#12, #7, #22, MT91 516, and CPR2002 Al 1) that are sufficiently 
bright and clearly separated from other bright sources for their 
spectra to be studied during the individual pointings. Table [3] 
summarizes the general properties of those X-ray sources that 
are associated with known OB stars. 

Figure[3] displays the hardness ratios HR2 = (H-M)/(H+M) 
versus HR1 = (M-S)/(M+S) for the OB-stars, where S, M, 
and H stand for the EPIC-pn count rates in the soft, medium, 
and hard energy bands, respectively. The solid lines yield the 
expected hardness ratios for single temperature plasma emis- 
sion at kT = 0.4, 0.6, and 0.85 keV for neutral hydrogen col- 
umn densities between 0.5 and 1.5 x 10 22 cirT 2 . The dashed line 
yields the average O-star X-ray spectrum as inferred by Naze 
(2009) and absorbed by column densities ranging between 0.5 
and 2.0 x 10 22 cm~ 2 . This latter model features two thermal 
plasma components with kT\ and kT2 equal to 0.24 and 2.3 keV, 
respectively, and the hotter component contributing 16% of the 
emission at 1 keV. As can be seen from this figure, with a few ex- 
ceptions, the hardness ratios of the majority of the O-type stars 
in Cyg OB2 can be understood in terms of typical O-star emis- 

4 At this stage, we stress that the star MT91 267 listed by Massey & 
Thompson ( TlWt and object CPR2002 All of Com eron et al. (2002) 
refer apparently to the same star. Kiminki et al. (2007) classified 
MT91 267 as a G-type star, whilst Negueruela et al. (2008) quoted an 
07.5 Ibf spectral type for All. Re-examination of the Kiminki et al. 
( 120071 ) Keck spectrum of MT91 267 by Kobulnicky (201 1, private com- 
munication) confirmed that this star is indeed an O-type star. 
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Fig. 3. EPIC-pn hardness ratios of the OB stars in Cyg OB2. The 
solid lines yield the synthetic hardness ratios of absorbed single 
temperature thermal plasma models for temperatures of kT = 
0.4, 0.6, and 0.85 keV (from bottom to top). The neutral column 
density increases along these lines from 0.5 to 1.5 x 10 22 cnT 2 
going from left to right. The dashed curve in turn indicates the 
hardness ratios of the 'typical' O-star two-temperature X-ray 
spectral model determined by Naze (2009) with an absorption 
column-density increasing from 0.5 to 2.0 x 10 22 cirT 2 . 



sion spectra affected by a large interstellar absorption. There is 
no evidence that the majority of the O-type stars have intrin- 
sically hard spectra. The hardest X-ray emission is observed 
for the B-type stars. However, it remains unclear whether these 
B-stars (except for CygOB2 #12) are intrinsic X-ray emitters. 
Their emission might also originate from hidden low-mass pre- 
main-sequence companions (e.g. Evans et al. 1201 11 and refer- 
ences therein), although we stress that in our dataset the B stars 
display little evidence of X-ray variability that would hint at the 
presence of a low-mass companion. 

From the results of the cross-correlation with the 2MASS 
data, we compiled near-IR colour-magnitude and colour-colour 
diagrams. For this purpose, we used the colour transformations 
of Carpenter (1200 11 updated in 2003^} to convert the 2MASS 
colours into the homogenized JHK photometric system of 
Bessell & Brett ( 11988I ). Figure|4] illustrates the results. The ma- 
jority of the massive stars that are X-ray sources (except for 
CygOB2 #12) have near-IR colours and magnitudes that are 
consistent with a reddening of Ak ~ 0.5 - 0.8 (see left panel 
of Fig. [4]). The locus of the fainter counterparts in this diagram 
is considerably shifted to the red, which implies that these stars 
are on pre-main-sequence (PMS) tracks. Bica et al. (2003]) re- 
ported average values of E(H - K s ) = 0.31 and A^ = 0.78 
for the two clusters they identified at the core of Cyg OB2. The 
majority of the near-IR counterparts of our X-ray sources are 
consistent with this value of reddening and suggest an age of 
2-3 Myr when compared to the PMS isochrones of Siess et al. 
(2000). These results are in excellent agreement with the conclu- 
sions of Albacete-Colombo et al. (2007a) and Wright & Drake 
(120091 ). although Wright et al. (120101) inferred a slightly older 

5 http://www.ipac.caltech.edu/2mass/index.html 
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Table 2. Properties of the X-ray sources in Cyg OB2 detected with XMM-Newton. 



[1] 


[2] 




[3] 


[4] 




[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


6 


20:32:22.3 


+41 


18:17.3 


0.1584D+00 


0.0162 


6 


0.6659D+00 


0.0000 


6 


0.2717D+00 


0.0000 


6 


13 


20:32:31.4 


+41 


14:06.0 


0.6859D- 


-02 


0.0193 


6 


0.5691D-01 


0.0000 


6 


0.4197D-01 


0.0000 


6 


16 


20:32:36.5 


+41 


22:13.1 


0.3771D- 


-03 


0.1853 


4 


0.2176D-02 


0.0000 


6 


0.1183D-02 


0.0000 


5 


22 


20:32:40.8 


+41 


14:28.0 


0.3500D- 


-01 


0.0000 


6 


0.4164D+00 


0.0000 


6 


0.3087D+00 


0.0000 


6 


27 


20:32:44.2 


+41 


24:22.3 


0.3257D- 


-03 


0.4902 


4 


0.1518D-02 


0.0083 


6 


0.8545D-03 


0.0017 


5 


34 


20:32:46.6 


+41 


18:06.8 


0.2761D- 


-03 


0.0000 


3 


0.8789D-03 


0.0000 


3 


0.6308D-03 


0.0000 


4 


74 


20:33:08.6 


+41 


13:17.8 


0.9955D- 


-02 


0.2646 


5 


0.5333D-01 


0.0707 


5 


0.1364D-01 


0.1639 


5 


80 


20:33:10.6 


+41 


15:07.6 


0.4022D- 


-01 


0.1605 


6 


0.2614D+00 


0.0005 


6 


0.1775D+00 


0.1212 


6 


92 


20:33:13.9 


+41 


13:04.4 


0.9154D- 


-03 


0.6731 


6 


0.4418D-02 


0.0000 


6 


0.8938D-03 


0.0000 


6 


95 


20:33:13.9 


+41 


20:21.1 


0.1232D- 


-01 


0.0157 


6 


0.4790D-01 


0.0814 


6 


0.9910D-02 


0.1750 


6 


96 


20:33:14.9 


+41 


18:50.0 


0.3530D+00 


1.0000 


6 


0.1292D+01 


0.0000 


6 


0.5778D+00 


0.0000 


6 


122 


20:33:23.5 


+41 


09:11.5 


0.9762D- 


-02 


0.1457 


6 


0.7324D-01 


0.0000 


6 


0.5054D-01 


0.0457 


6 


134 


20:33:26.6 


+41 


10:58.4 


0.8776D- 


-03 


0.0082 


5 


0.543 1D-02 


0.0000 


6 


0.2440D-02 


0.0000 


6 



[1] 


[13] 




[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 




[23] 


[24] 


6 


0.4100D- 


-01 


0.0001 


6 


0.2546D+00 


0.0000 


6 


0.1081D+00 


0.0000 


6 


0.4230D- 


-01 


0.0020 


6 


13 


0.1747D- 


-02 


0.2114 


6 


0.2484D-01 


0.0000 


6 


0.1794D-01 


0.0000 


6 


0.1948D- 


-02 


0.0024 


6 


16 


0.1344D- 


-03 


0.6146 


4 


0.7055D-03 


0.0000 


6 


0.6065D-03 


0.0000 


5 


0.8757D- 


-04 


0.6374 


3 


22 


0.8892D- 


-02 


0.0032 


6 


0.1429D+00 


0.0000 


6 


0.1087D+00 


0.0000 


6 


0.7893D- 


-02 


0.0000 


6 


27 


0.1216D- 


-03 


0.4100 


2 


0.5925D-03 


0.3433 


4 


0.4596D-03 


0.9076 


3 


0.1947D- 


-03 


0.9668 


4 


34 


0.1740D- 


-03 


0.3820 


3 


0.7237D-02 


0.0000 


2 


0.3863D-03 


0.0000 


4 


0.1304D- 


-03 


0.0282 


2 


74 


0.1779D- 


-02 


0.0166 


6 


0.1772D-01 


0.1740 


6 


0.5399D-02 


0.0019 


6 


0.1847D- 


-02 


0.1020 


6 


80 


0.9531D- 


-02 


0.0093 


6 


0.8768D-01 


0.0370 


6 


0.6220D-01 


0.0007 


6 


0.9659D- 


-02 


0.0000 


6 


92 


0.2085D- 


-03 


0.6109 


5 


0.7190D-03 


0.0879 


6 


0.4920D-03 


0.0000 


5 


0.2043D- 


-03 


0.8327 


5 


95 


0.2614D- 


-02 


0.0900 


6 


0.1552D-01 


0.2145 


6 


0.3727D-02 


0.7770 


6 


0.2804D- 


-02 


0.1776 


6 


96 


0.8258D- 


-01 


1.0000 


6 


0.4440D+00 


0.0000 


6 


0.2159D+00 


0.0000 


6 


0.8393D- 


-01 


1.0000 


6 


122 


0.244 1D- 


-02 


0.2976 


6 


0.2709D-01 


0.0004 


6 


0.1681D-01 


0.0000 


6 


0.2605D- 


-02 


0.9902 


6 


134 


0.2774D- 


-03 


0.4752 


6 


0.1713D-02 


0.0034 


6 


0.6320D-03 


0.0000 


5 


0.2919D- 


-03 


0.7785 


4 



[1] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


[35] 


[36] 


[37] 


6 


0.2758D+00 


0.0000 


6 


0.1158D+00 


0.0000 


6 








5.187 


4.745 


4.339 


AEA 


13 


0.2564D-01 


0.0000 


6 


0.1902D-01 


0.0000 


6 








7.817 


7.094 


6.664 


AAA 


16 


0.4226D-03 


0.0000 


6 


0.6045D-03 


0.0000 


2 








12.106 


11.337 


10.879 


AAA 


22 


0.1431D+00 


0.0000 


6 


0.1080D+00 


0.0000 


6 




20.104 


19 


4.667 


3.512 


2.704 


DDD 


27 


0.5709D-03 


0.6657 


5 


0.1174D-02 


0.0081 


5 








13.820 


12.754 


12.408 


AAA 


34 


0.4953D-03 


0.0000 


4 


0.465 1D-03 


0.0000 


6 








16.337 


14.868 


14.474 


AAA 


74 


0.1710D-01 


0.0000 


6 


0.5112D-02 


0.0005 


6 


20.747 


16.699 


95 


7.110 


6.540 


6.226 


AAA 


80 


0.9184D-01 


0.0000 


6 


0.6456D-01 


0.0004 


6 


19.642 


15.811 


103 


6.468 


5.897 


5.570 


AAA 


92 


0.1364D-02 


0.0001 


6 


0.1356D-02 


0.0000 


5 


20.904 


17.445 


118 


9.034 


8.559 


8.280 


AAA 


95 


0.1608D-01 


0.2545 


6 


0.3704D-02 


0.5000 


6 


17.385 


14.174 


120 


7.248 


6.818 


6.611 


AAA 


96 


0.4518D+00 


0.0000 


6 


0.2164D+00 


0.0000 


6 


15.267 


12.462 


126 


6.123 


5.721 


5.503 


AAA 


122 


0.2510D-01 


0.0965 


6 


0.1775D-01 


0.3997 


6 








7.025 


6.380 


6.050 


AAA 


134 


0.1488D-02 


0.1372 


6 


0.9622D-03 


0.0000 


5 








8.971 


8.434 


8.165 


AAA 



[1] 


[38] 


[39] 


[40] 


[41] 


[42] 


6 








06.5-71 + OB-Ofpe/WN9 


CygOB2#5 


13 


12.87 


1.22 


2.19 


07.5M 


CPR2002 All 


16 


16.19 


1.16 


2.06 




MT91 291 


22 


11.46 


2.37 


3.35 


B5Ie 


CygOB2#12 


27 












34 












74 


11.55 


0.74 


2.04 


04111(f) 


CygOB2#22 


80 


10.96 


0.69 


1.81 


05 If + 06-7 


CygOB2#9 


92 


12.92 


0.54 


1.81 


08 V 


MT91 455 


95 


10.55 


0.19 


1.45 


03 If 


CygOB2#7 


96 


9.06 


0.08 


1.30 


05.51(f) + 06 


CygOB2#8a 


122 


11.84 


0.93 


2.20 


05.5 V 


MT91 516 


134 


13.00 


0.61 


1.87 


08.5V 


MT91 534 



Notes. This table will be made available in its entirety at CDS. Some excerpts are shown here for guidance. The columns have the following 
meaning: [1] = source number; [2] and [3] = source coordinates (J2000.0); [4] - [6] = mean count rate in EPIC-pn soft band, probability that the 
source be constant and number of observations of this source in this band with EPIC-pn; [7] - [9] = same for the EPIC-pn medium band; [10] - 
[12] = same for the EPIC-pn hard band; [13] - [15] = same for the EPIC-MOS1 soft band; [16] - [18] = same for the EPIC-MOS1 medium band; 
[19] - [21] = same for the EPIC-MOS1 hard band; [22] - [24] = same for the EPIC-MOS2 soft band; [25] - [27] = same for the EPIC-MOS2 
medium band; [28] - [30] = same for the EPIC-MOS2 hard band; [31] and [32] = UVM2 and UVWl magnitudes; [33] = ID number of OM 
counterpart; [34] - [36] = J, H, K s photometry from 2MASS; [37] = 2MASS quality flag or "ncp", "2cp" "3cp" where 2cp (3cp) means there are 
2 (resp. 3) 2MASS counterparts inside the 4 arcsec correlation radius whilst ncp means no 2MASS counterpart was found inside the correlation 
radius; [38] - [40] = V, U — B, B — V from Massey & Thompson ( 1991 1; [41] and [42] = spectral type and common identifier. 
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Table 3. General properties of the OB stars in the core of Cyg OB2 as seen with XMM-Newton. 



Name 


Spectral type 






pn 




MOS1 


MOS2 


UVM2 


UVWl 






CR 


HR1 


HR2 


var 


CR 


var 


CR 


var 






Objects studied in previous papers 

Cyg OB2 #5 06.5-7 1 + OB-Ofpe/WN9 


1.096 


0.62 


-0.42 


SMH 


n A C\A 

U.4U4 


OA /fT_T 

.SMrt 


n a i a 
0.434 


aJVllrl 






Cyg OB2 #8a 


05.51(f) + 06 


Z.ZZJ 


U.J / 


U.jo 


ivlrl 


0.742 


MH 


0.752 


MH 


1 J.Z / 


1 Af* 
lZ.<+0 


CygOB2#9 


05 If + 06-7 


u.4/y 


n 7^ 

U. / J 


n 1 q 


1V1 


0.159 


SMH 


0.166 


SMH 


1 Q f*A 


1 ^ C 1 


X-ray bright OB stars analysed in this paper 
Cyg OB2 #7 03 If 


0.070 


0.59 


-0.66 


s 


0.022 




0.023 




17.39 


14.17 


Cyg OB2#12 


B5Ie 


0.760 


0.84 


-0.15 


SMH 


0.260 


SMH 


0.259 


SMH 




20.10 


CygOB2#22 


04111(f) 


0.077 


0.69 


-0.59 




0.025 


SH 


0.024 


MH 


20.75 


16.70 


MT91516 


05.5 V 


0.133 


0.76 


-0.18 


MH 


0.046 


MH 


0.045 








CPR2002A11 


07.5 Ibf 


0.106 


0.78 


-0.15 


SMH 


0.044 


MH 


0.047 


SMH 






X-ray fainter OB stars, not analysed in detail 






















Cyg OB2 #4 


07 Illf 


0.018 


0.23 


-0.25 




0.011 


H 


0.010 








Cyg OB2 #6 


08 V 


0.009 


0.04 


-0.89 




0.002 




0.002 




16.79 


13.82 


Cyg OB2 #8c 


04.5 m(fc) 


0.040 


0.42 


-0.78 


S 


0.014 




0.014 




16.58 


13.41 


CygOB2#15 


08 V 


0.009 


0.01 


-0.63 


MH 


0.003 


S 


0.003 








Cyg OB2#16 


07 V 


0.009 


0.28 


-0.79 




0.003 




0.003 




16.72 


13.90 


Cyg OB2#17 


08 V 


0.005 


0.09 


-0.57 




0.002 




0.002 




18.18 


15.13 


Cyg OB2 #24 


07 V 


0.009 


0.41 


-0.59 


S 


0.002 




0.003 




19.39 


15.97 


Cyg OB2 #27 


09.5 V 


0.003 


0.26 


-0.47 


H 


0.001 




0.001 




19.82 


16.45 


MT91 311 


B2V 


0.003 


0.61 


0.28 




0.001 




0.001 




20.94 


18.01 


MT91 322 


B2.5V 


0.002 


0.48 


-0.13 




0.001 


H 


0.001 






18.71 


MT91 336 


B3III 


0.004 


0.73 


-0.18 




0.001 




0.001 




20.55 


17.75 


MT91376 


08 V 


0.005 


0.12 


-0.59 




0.001 




0.001 




18.44 


15.43 


MT91448 


06 V 


0.004 


0.55 


-0.89 




0.002 




0.001 


M 




19.04 


MT91455 


08 V 


0.006 


0.66 


-0.66 


MH 


0.001 


H 


0.003 


MH 


20.90 


17.45 


MT91485 


08 V 


0.007 


0.43 


-0.65 




0.002 




0.003 




19.25 


16.03 


MT91534 


08.5 V 


0.009 


0.72 


-0.38 


SMH 


0.003 


MH 


0.003 


H 






MT91573 


B3I 


0.002 


0.51 


-0.19 


MH 






0.001 






18.12 



Notes. Columns 3, 7, and 9 provide the average count rates (cts s -1 ) over the 0.5 - 8.0 keV energy band for the pn, MOS1, and MOS2 detectors, 
respectively. Columns 6, 8, and 10 indicate the energy bands for which the null hypothesis of a constant count rate is rejected at the 5% level for the 
pn, MOS1, and MOS2 detectors. The hardness ratios (quoted in columns 4 and 5 for the pn instrument only) are defined as HR1 = (M-S)/(M+S) 
and HR2 = (H-M)/(H+M). The X-ray spectra of the first three objects were discussed in previous papers, those of the five following objects are 
discussed in detail in this paper and the remainder of the objects are too faint for spectral analysis. 




Fig. 4. Left panel: (J - K,K) near-IR colour-magnitude diagram of the 2MASS counterparts of our X-ray sources. The main- 
sequence (from Knodlseder 2000) is shown for a distance modulus of 10.4 and for two different values of the reddening (from left 
to right, Ak = 0.0 and 0.78, Bica et al. 2003). The reddened pre-main sequence isochrones from Siess et al. (120001 1 are shown for 
ages of 2 and 3 Myrs (dashed lines). Right panel: (H - K, J - H) colour-colour diagram of the 2MASS counterparts. The solid lines 
yield the locus of the main-sequence according to Bessell & Brett ( 1988), while the dashed straight line indicates the locus of the 
unreddened classical T-Tauri stars following Meyer et al. ( 1997 ). The reddening vector is shown according to Bica et al. (2003 ). 
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Fig. 5. Left panel: colour-magnitude diagram compiled from the OM photometry of CygOB2. Right panel: (J - K,K) near-IR 
colour-magnitude diagram (as Fig.|4ji of the 2MASS counterparts of the OM sources that are detected in both UV filters. The F-type 
star indicated in these figures is Cyg OB2 #72 (= MT91 664, Massey & Thompson \l 99 II) . 



age of 3.5^jQgMyr. We draw attention to a group of sources 
that fall roughly along the main-sequence relation (shifted to the 
distance of Cyg OB2 and reddened according to the Bica et al. 
2003 results), but below the main-sequence turn-on. One of these 
sources (number 86) displays a decaying flare during observa- 
tion 1. These objects are likely to be late-type active foreground 
stars unrelated to CygOB2. From the right panel of Fig.|4] we 
see that there is little evidence of large IR excesses as one would 
expect for classical TTauri stars (see Meyer et al. 119971 . This is 
in line with the conclusion of Albacete-Colombo et al. (f2007a ) 
and Wright et al. ( 120101 1, who found that only 4 - 8% of the near- 
IR counterparts of their X-ray sources display infrared excesses 
that could hint at the presence of an accretion disk. 

Vink et al. (2008) performed a search for Ha emitters in 
Cyg OB2. Their field C overlaps to a large extent with our EPIC 
field of view. Four Ha emitters tabulated by Vink et al. (2008 ) lie 
inside the EPIC FoV. Only one of them is marginally consistent 
with a detection in our X-ray data (source C9]Q although we do 
not regard this correlation as a significant detection. 

Using our OM data, we also performed a search for UV 
sources. After rejecting some spurious sources (mainly due to 
the artifacts mentioned in Sect. |2]i, we found a list of 254 objects 
that are detected at least in the UVWl filter and at least during 
two different observations. Using a cross-correlation radius of 
4 arcsec, we find that 32 X-ray sources have a counterpart in our 
list of OM sources. Out of these 32 objects, 19 are measured in 
both OM filters, whilst the others are too faint in the UVM2 filter 
to be detected. Among the X-ray sources detected in both UV fil- 
ters, 15 have an OB-type counterpart. Six X-ray flaring sources 
fall inside the OM FoV. Only one of them (the X-ray source 139) 
is detected with a magnitude of 18.8 in UVWl (no detection in 
UVM2). From a cross-correlation with the catalogue of Massey 
& Thompson (1991), we find counterparts for 225 of our OM 
sources within a correlation radius of 2.0 arcsec. Adopting the 
same correlation radius, all but five OM sources are found to 

6 We note here that the position of source C9 in Fig. 1 of Vink et al. 
(2008) is inconsistent with the coordinates given in their Table 2. We 
used the latter in our cross-correlation. 



have a 2MASS counterpart. The results of this correlation, along 
with the results of a x 2 variability test, are provided in Table 4 
that is available at the CDS. 

The UVM2 filter contains the broad interstellar absorption 
feature at 2200 A and is thus quite sensitive to the interstellar col- 
umn density (Royer et al. 120001 ). Given the rather heavy redden- 
ing towards the Cyg OB2 association, objects that are detected in 
both filters correspond to either the brightest and hottest associ- 
ation members or foreground stars. In the UV colour-magnitude 
diagram (see Fig. [5]), the O-type stars define a clear, almost ver- 
tical sequence around UVM2 - UVWl - 3. The latter value is 
consistent with the large reddening towards the Cyg OB2 asso- 
ciation (Royer et al. 120001 ). 

For the O-type stars in our OM FoV, we inspected the inter- 
pointing light curves of these stars that were found to be vari- 
able using a x 2 test whilst not exceeding the limiting brightness 
at which 10% coincidence losses occurr. Among the stars that 
did not display significant UV magnitude variations, the most 
remarkable cases are CygOB2 #7, 9, and 12. The O-type stars 
that were found to be variable in either of the UVM2 and U VW1 
bands are, respectively, CygOB2 #8a, 8c, 16, 27 and CygOB2 
#22, 23, 27, MT91421, and 507. However, in most cases the 
peak-to-peak amplitude of this variability is limited to less than 
0. 15 mag. Whilst this is larger than the formal errors in the mag- 
nitudes, the problems caused by the straylight features probably 
account for part of this apparent variability. There is one excep- 
tion though: Cyg OB2 #27 (Table|5]). The latter is a W UMa con- 
tact binary consisting of a late O and an early B-type star in a 
1.46 day orbit (Rios & DeGioia-Eastwood 2004 ). Unfortunately, 
there are no ephemerides available for this system, but we sus- 
pect that the 0.64 mag drop that we observe in the UVM2 filter 
is likely due to the eclipse of the O-type star by its cooler com- 
panion. 

3.2. Long and short term variations 

Using a x 1 test, we searched for inter-pointing variability. The 
variability test is performed using the mean exposure-corrected 
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Table 5. OM light curve of Cyg OB2 #27 (= MT91 696). 



JD-2450000 


AB-magnitude 




UVW1 






3308.5052 


±0.0184 


16.38 


±0.01 


3318.4781 


± 0.0242 


16.41 


±0.01 


3328.4983 


± 0.0289 


16.55 


±0.01 




UVM2 






3308.6205 


± 0.0463 


19.86 


±0.11 


3318.6049 


±0.0521 


19.63 


±0.08 


3328.5934 


± 0.0579 


20.27 


±0.13 


3338.5573 


± 0.0473 


19.78 


±0.08 



Notes. The quoted uncertainty in the time corresponds to half of the 
exposure time. 




1 7.5 




20.5 r 

-0.5 0.5 1 1.5 

Phase (P = 1.46 days) 

Fig. 6. OM light curves of Cyg OB2 #27 as a function of orbital 
phase. Phase 0.0 was arbitrarily chosen at the time of the first 
observation. 



EPIC count rates from the different pointings, along with their 
estimated errors. The results are summarized in Table|6]for an 
adopted significance level of 95%. As can be seen from that ta- 
ble, the incidence of variability increases significantly from the 
soft to the hard energy band. There are two reasons for this trend: 
(1) the huge interstellar absorption leads to rather low count rates 
in the soft energy band; and (2) a fraction of the variable sources 
display flares (see below) associated with a hardening of their 
spectra, hence the largest variability occurs in the hard energy 
band. 



Table 6. Statistics of our inter-pointing x 2 variability test as a 
function of instrument and energy band. 



Instrument 


Soft 


Medium 


Hard 


Full 


MOS1 


5 


29 


39 


52 


MOS2 


5 


27 


39 


48 


pn 


14 


43 


47 


69 



Notes. The numbers provide the fraction (in %) of sources that display 
variability at a significance level of > 95%. 



The variability comes in different flavours. In some cases, 
it is clearly associated with a large (a factor 10 or more) in- 
crease in the count rate during a specific observation (sources 
number 34, 39, 67, 132, and 173). In other cases, the count 
rates display a modulation of typically a factor 2-3 in certain en- 
ergy bands (sources number 59, 107, 126, 144, 160, 182) from 
one observation to the other. In yet another group of (bright) 
objects, the variability is even more subtle, but still significant 
(CPR2002 All). We note that the list of variable sources, espe- 
cially those of the last category, is not restricted to faint (proba- 
bly low-mass) sources, but also includes OB-type stars. Table|3] 
indicates that 9 OB stars out of 25 were detected to be vari- 
able, at least in one energy band, by at least two of the three 
EPIC instruments. This is in line with the conclusion of Naze 
(2009), who found that between one-third and two-thirds of the 
massive stars in the 2XMM catalogue display long-term varia- 
tions. In our sample of massive stars, the variable sources in- 
clude known colliding-wind binary systems (Cyg OB2 #5, #8a, 
#9), radial velocity variables (MT91 516, Kiminki et al. l2"00"7T l. 
as well as stars for which multiplicity remains currently undeter- 
mined (Cyg OB2 #12, CPR2002 All, Cyg OB2 #22, MT91 455, 
MT91 534). Among the eight X-ray brightest O-stars (those with 
mean pn count rates above 0.05 cts s~'), the only object that does 
not display significant variations in our data is Cyg OB2 #7. This 
star in addition does not display significant radial velocity vari- 
ations (Kiminki et al. 120071) . These results suggest that in our 
sample the X-ray variability of O-type stars could be an indica- 
tor of multiplicity. 

Adopting a timebin of 1 ksec, we extracted the light curves 
of those sources that display significantly higher count rates dur- 
ing a specific observation for each instrument and over three dif- 
ferent energy bands (medium, hard, and total). A priori, these 
long-term variations could represent a flaring event during the 
observation with the highest count rate. X-ray flaring is a com- 
mon property of magnetically active stars (e.g. Favata et al. 
120051) . Albacete Colombo et al. d2007bl l detected flaring activity 
in about 15% of the sources in their 100 ksec Chandra observa- 
tion of Cyg OB2. 

The X-ray light curves of the flare-candidates actually dis- 
play a wide variety of behaviours (see Fig.[7]i: 

• Eight sources (34, 39, 64, 67, 86, 132, 149, and 173) display 
a genuine flare, although in two cases (64 and 86), our ob- 
servations only cover the decay phase of an event that most 
likely started before our observation. We attempted to fit the 
decay phase of the light curves in the total energy band with 
an exponential decrease of the form CR max exp [— ij, where 
CR m ax i s tr, e count rate at the maximum of the flare and r is 
the exponential decay time. The resulting decay times of five 
objects are between 1.7 and 7.3 ksec (i.e. 0.5 and 2 hours). 
Albacete-Colombo et al. (2007b) found decay times rang- 
ing from 0.5 to 10 hours in their 100 ksec observation. Our 
dataset is clearly biased against the detection of flare events 
with decay times longer than about 25 ksec (i.e. the mean 
duration of our individual observations). We note that some 
flare events cannot be modelled by an impulsive rise fol- 
lowed by an exponential decay. Source 173 for instance rises 
rather slowly (within 10 ksec) to its maximum and remains 
there for the remainder of observation 5. The flare event of 
source 34 is probably an example of sustained heating with 
a second flare occuring about 6 ksec after the first, more en- 
ergetic, one (see Fig. [8]). 

• Twenty-five sources (27, 57, 68, 70, 78, 82, 84, 89, 93, 94, 
101, 106, 109, 124, 128, 130, 135, 137, 139, 140, 153, 179, 
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Fig. 7. EPIC light curves of variable sources illustrating the variety of behaviours that are observed. The green, blue, and black 
histograms yield the background-subtracted light curves in the medium, hard, and total energy band respectively. The time is given 
in ksec from the start of the corresponding observation. Top, from left to right: sources 39, 132, and 190 during observations 4, 6, 
and 2, respectively. Bottom, from left to right: sources 67 and 176 during observations 1 and 2, respectively. 



186, 190, and 196) display a light curve with no clear flare 
events, but with some short-term (a few ksec) variations. In 
most cases, a significant fraction (typically more than half) 
of the photons of these sources are emitted in the hard energy 
band. Two of these sources (70 and 78) have a higher count 
rate during two observations (3 + 5 and 2 + 5, respectively). 
• Eight sources (16, 73, 79, 92, 134, 164, 169, and 176) dis- 
play a high count rate during a specific observation with slow 
variations (i.e. on the timescale of the duration of the obser- 
vation). These objects could be candidates for long duration 
(significantly longer than the duration of the observation) 
flares. 

3.3. Spectra of variable non-OB sources 



As pointed out in Sect. 3.2 a number of sources appear much 



brighter during a single observation than in average, although 
not all of these sources actually display a genuine flaring be- 
haviour. For the brightest of these objects, we were able to ex- 
tract and analyse their EPIC spectra. We attempted to fit these 
spectra with first an absorbed thermal plasma model and then an 
absorbed power-law model using the xspec software (version 
12.6.0, Arnaud fl9961 l. Except for a few cases (sources 73, 176) 
where both models provide equally good fits, the absorbed ther- 
mal plasma model (wabs x apec) is the one that most closely fits 
the data. The results are summarized in Tableland the fit of the 
spectra of source 34 is illustrated in Fig. [8] 

All these sources have rather hard spectra characterized by 
high plasma temperatures (kT > 4keV) as expected for mag- 
netically triggered X-ray flares. The column densities are some- 
what lower than expected for Cyg OB2 members, especially in 
the case of source 53^] However, we stress that, given the high 

7 We do not have a 2MASS counterpart with reliable photometry for 
this source, but both the OM photometry and the optical photometry 



plasma temperatures, the column densities are rather poorly con- 
strained. We emphasize that the quoted fluxes correspond to av- 
erage values over the entire observation under consideration. For 
flaring sources, they do not represent the peak luminosity. We 
note that the errors in the fluxes should be regarded only as in- 
dicative since the relevant xspec routine does not account for 
possible correlations between model parameters. 



4. X-ray spectra of OB-type stars 

At the sensitivity level of the current generation of X-ray tele- 
scopes, most early-type stars do not display short-term X-ray 
variability. However, in massive binary systems and single mas- 
sive stars with magnetic fields of sufficient strength to confine the 
wind in the equatorial region, a modulation of the X-ray emis- 
sion is expected and has been observed on the orbital or rota- 
tional period typically of the order of either days or weeks (see 
Giidel & Naze 2009 and references therein). 

The X-ray spectra of three of the brightest sources, 
(Cyg OB2 #5, #8a, and #9) were discussed elsewhere (De Becker 
et al.|2006l Linder et al. |2009l Naze et al. |2010l and Blomme et 
al. 1201 Oi l, and we do not repeat their analysis here. For the other 
OB-stars that are sufficiently X-ray bright, we analysed the spec- 
tra of each observation using xspec. 

The interstellar neutral-hydrogen column density is evalu- 
ated from the E(B - V) colour excess computed from the B - V 
colour measured by Massey & Thompson (1991) and by adopt- 
ing the gas-to-dust ratio of Bohlin et al. d!9781 l. Early-type stars 
host dense stellar winds that can contribute to the absorption of 
X-rays produced in the inner regions of these winds. Therefore, 
whenever this turned out to be supported by our observations, we 



from MT91 support the idea that this star is less heavily reddened than 
the majority of the stars towards Cyg OB2. 
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G. Rauw: A multi-epoch XMM-Newton campaign on Cyg OB2 
Table 7. Results of the spectral fits of the flaring sources with single temperature thermal plasma models. 



wabs x apec(T) 



Object 


Obs. 


Inst. 


d.o.f. 


r 


N H 
HO 22 crrT 2 1 


kT 

flceVI 
\AC v ) 


norm 

HO -14 rm- 5 ) 


fx 

( pro nm~ 2 s — ^ 




fx 

i -2 
i pro pm - 


) 


67 


1 


M1+M2 


21 


0.92 


1 22 + - 75 


12 2 


(3.6 +, J) x 10- 4 


(5.24 + Hp x 10- 


13 


7.32 x 10" 


13 


73 
176 


2 
2 


EPIC 
M1+M2 


11 

38 


1.27 
1.22 


0.90!;| 


A 7+6-3 
-L6 

6 2 +12 - 9 

-2.1 


(1.0 +4 )x IO" 4 
(11.0ff|)x io- 4 


(1.06+™) x 10- 
(13.83 + ™ 4 ) x 10 


13 
-13 


1.64 x 10- 
20.46 x 10 


13 

-13 


34 


3 


EPIC 


87 


0.78 


1.02 +: | 2 


8.ir? : | 


(4.4 + 4 ) x 10~ 4 


(6.17 + !!- 22 )x 10- 


13 


8.68 x 10- 


13 


53 


3 


pn 


30 


0.75 


26 + ( " 


5.5! 2 9 3 


(1.9 +:2 ) x 10- 4 


(2.86 + ° :24 ) x 10- 


13 


3.43 x 10" 


13 


39 
173 


4 
5 


EPIC 
EPIC 


34 
29 


1.15 
1.49 


1.17 + -g 
1.1 3 +:3 ' 


9.6;^ 
6 9 +17 - 8 


(3.1 +: ^)x io- 4 
(6.1 + j^)xl0- 4 


x 10- 
(7.97 + "p x 10" 


13 
13 


6.16 x 10- 
11.62 x 10 


13 
-13 


132 


6 


Ml +M2 


25 


0.76 


1.08+f 




(4.5 + £°) x IO" 4 


(5.94 + ° :4 [|) x 10- 


13 


8.62 x 10- 


13 



Notes. The last two columns yield the observed and absorption corrected flux over the 0.5 - 10.0 keV energy range. 
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Fig. 8. Left: EPIC spectra (along with the best-fit model) of source 34 during the third observation. Right: EPIC light curves of this 
source during the same observation. The green, blue, and black histograms indicate the light curves in the medium, hard, and total 
energy bands, respectively. The time is given in ksec elapsed from the start of the third observation. 



included an ionized wind absorption model (Naze et al. 2004]> in 
addition to the neutral interstellar-medium absorption. 

We compare the properties of these OB stars to the general 
X-ray emission properties derived by Naze 020091 1, who anal- 
ysed the data of OB stars in the 2XMM catalogue. Naze ((2009) 
obtained log ^ = -6.45 ± 0.51 for those O-stars that have data 

of sufficient quality for a spectral fit and log j^- = -6.97 ± 0.40 

Mm] 

from 26 O-type stars in the region of Cyg OB^j In cases where 
a two-temperature model was required to fit the 2XMM EPIC 
spectra, the normalization of the hotter component was gener- 
ally found to be less than 0.3 times that of the lower temperature 
component (Naze |2009l . 



8 Note that the latter relation was obtained by combining the data 
for stars for which the X-ray fluxes were either extracted from spectral 
fits or obtained from the EPIC count rates using an average conversion 
factor. Note also that these relations include both single and binary 
O-stars as Naze ( 2009 1 found no clear evidence of a general X-ray over- 
luminosity of O + O binary systems. Some specific interacting-wind O- 
star binaries are indeed X-ray overluminous, but this overluminosity is 
not a general feature and is often restricted to specific orbital phases. 



4.1. CygOB2#12 

Massey & Thompson dl991l) quote a colour of B - V = 3.35 for 
CygOB2 #12, which corresponds to a colour excess of E(B - 
V) - 3.45, assuming a B5 1 spectral type (see below). When we 
adopt a distance modulus of 10.4 mag and Ry = 3.1, we find 
that M v = -9.64. Assuming a bolometric correction of -0.95 
then leads to log ^ = 6.13. This luminosity places the star 
very near the empirical Humphreys-Davidson limit (Humphreys 
& Davidson [T984l . 

The reddening of CygOB2 #12 is larger than the average 
value of other OB stars in CygOB2, and the membership of 
this star is still uncertain. ESA's forthcoming GAIA mission 
should provide accurate astrometric measurements for this star 
that will help establish its membership. Spectral classifications 
of Cyg OB2 #12 range from B3 la to B8 la, suggesting that the 
star actually changes its spectral type with time (Kiminki et al. 
I20071 and references therein). For instance, Kiminki et al. (2007) 
presented spectra that indicate a B3 Iae type in September 2000 
and B6-8 one year later. The nature of this star remains contro- 
versial, although based on its likely spectroscopic variability and 
incredibly high luminosity, CygOB2 #12 has been suggested 
to be an LBV candidate (see Clark et al. 2005 and references 
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therein). Cyg OB2 #12 was detected as an IRAS source and is ap- 
parently surrounded by both warm and cold dust (Parthasarathy 
et al. [1^921 . 

Souza & Lutz ( fT980b classified CygOB2 #12 as a B8 su- 
pergiant with weak Balmer lines and emission in Ha suggesting 
that the other Balmer absorptions could be filled-in by emission. 
They also reported a variable absorption bluewards of Ha that 
they tentatively interpreted as a blueshifted Ha absorption as- 
sociated with an expanding shell at a velocity of 1400 km s~'. 
Klochkova & Chentsov d2U04l l classified Cyg OB2 #12 as B5 Ia + 
with an uncertainty of half a subtype. These authors presented 
evidence of a line radial-velocity gradient that might indicate an 
infall of matter. They estimated that the star has a radius of or- 
der 338 R and a mass-loss rate of 4 x 10~ 5 M yr~' . Klochkova 
& Chentsov (2004) inferred a rather slow wind velocity of 
150kms _1 , which would have important consequences for the 
understanding of the observed X-ray emission (see below). Such 
a slow wind would indeed be able to produce hard X-ray emis- 
sion neither by means of wind embedded shocks nor head-on 
collisions of the wind from two hemispheres in a magnetically 
confined wind model. If CygOB2 #12 were a binary system, 
then the wind of the secondary star, colliding with that of the B 
supergiant, would have to be much faster. We note however that 
Leitherer et al. (1982 ) derived a significantly larger wind veloc- 
ity of 1400km s _1 from the Ha line profile and that Wang & Zhu 
( 2003 ) reported an expanding shell at a velocity of 3 1 00 km s~ 1 . 

The EPIC spectra of this star cannot be fitted with 
single-temperature plasma models and we instead used two- 
temperature models. The huge interstellar column density that 
we infer from the E(B - V) colour excess (Nh.ism = 2.0 x 



10 cm ) produces such a strong absorption at low photon en- 
ergies that there is no room left for an additional ionized ab- 
sorption column associated with the stellar wind. After testing 
several options of absorbed two plasma temperature models, we 
found that the optimal results are obtained when only interstel- 
lar absorption is included and the absorbing column is left as 
a free parameter. The results of the fits are quoted in TableJH] 
and illustrated in Fig. [9] The average absorption-corrected X-ray 
flux corresponds to an X-ray luminosity of 4.1 x 10 33 ergs _1 . 
With the caveat that the large interstellar reddening renders the 
bolometric luminosity estimation and the X-ray flux correc- 
tion more uncertain than for other stars in our sample, we de- 
rive log y^- — -6.10. This is an unusually large ratio, signif- 
icantly larger than the average value for O-type stars and cer- 
tainly exceptional for B-type stars. These characteristics imply 
that Cyg OB2 #12 might either be a binary or host a stellar wind 
confined by a strong magnetic field. 

Figure[9] reveals moderate (~ 10%) flux variability on 
timescales from a few days to a few weeks, as well as large varia- 
tions (37%) on timescales of years^] Long-term variations in the 
column density are also seen. However, there are apparently little 
or no changes in the plasma temperatures which are found to be 
0.76 ± 0.03 and 2.03 + 0.19keV for the cooler and ho tter com- 
ponent respectively. Albacete-Colombo et al. (2007a) reported 
a roughly linear decay of the ACIS-I countrate of Cyg OB2 #12 



9 As a complement to our XMM-Newton data, we analysed the ASCA- 
SIS spectra of Cyg OB2 #12 from 29 April 1993 previously analysed by 
Kitamoto & Mukai (1996). The SIS spectra can be fitted with a single 
temperature component at kT = 0.90. The observed flux was found to 
be 2.2 x 10~' 2 ergcnr 2 s -1 , quite similar to the value of our observa- 
tion 6. However, owing to the modest angular resolution of ASCA and 
the contamination of the background by straylight from CygX-3, the 
uncertainty in this result is rather large (about 25%). 



from 0.1 8 to 0.16 cts s _1 (i.e. a 12% variation) over their lOOksec 
Chandra observation. This result is very much in line with the 
~ 18% variations in the observed flux that we measure in our 
first four data points on timescales of about a week and the 10% 
variations seen between the last two pointings. 



Table 9. OM UVWl data of Cyg OB2 #12. 



JD-2450000 


UVWl 


3308.505 


20.08 


±0.10 


3318.478 


20.13 


±0.09 


4220.282 


20.00 


±0.13 


4224.061 


20.25 


±0.16 



CygOB2 #12 was inside the OM FoV during four (1, 2, 5 
and 6) out of the six observations. Owing to its heavy reddening, 
the star is detected only in the UVWl band (see Table|9]) with an 
average magnitude of 20. 1 1 . The lack of significant variations 
in this magnitude with time, indicates that if changes in spectral 
type have occured during our campaign, they were not associated 
with large changes in the UVWl spectral domain. 

4.2. CygOB2#7 

This 03 If supergiant displays rather soft X-ray spectra that 
can be fitted well using a single temperature plasma compo- 
nent. From the colour excess E(B — V) = 1.73 (Massey & 
Thompson 119911 ). we infer an interstellar column density of 
A^ M = 1.0 x 10 22 cirT 2 . The spectra show evidence of additional 
absorption and we thus included absorption by an ionized wind 
in our model. The results of the fit (at energies above 0.4 keV) 
are listed in Table 10 There is no strong evidence of variability. 
The observed flux and the wind column density remain constant 
to within 10%. At first sight, the plasma temperature varies be- 
tween 0.56 and 0.77 keV. However, there is some degeneracy in 
the fits between the plasma temperature and its normalization 
and all the spectra are actually consistent with each other. 

Herrero et al. (120021 ) inferred a luminosity of log ^ =5.91 + 
0.10 (assuming the distance modulus of Massey & Thompson 
1 1 99 lb - Scaling this number to the lower distance modulus of 
Hanson (f2003, DM = 10.4) that we have adopted in this pa- 
per, yields log ^ = 5.59. Proceeding in the same way with the 
bolometric magnitude derived by Negueruela et al. (2008 ) yields 
log ^ = 5.87. Our mean, ISM absorption-corrected X-ray flux 

(8.76 + 1.26)10 



13, 



,-2 o-l 



'erg cm "s 1 yields log j^- = -6.99 + 0.06 
or -7.27 ± 0.06 for the Herrero et al. ( 120021 ) and Negueruela 
et al. ( 120081) bolometric luminosities, respectively. We conclude 
that the X-ray emission of Cyg OB2 #7 does not reveal any be- 
haviour (hard X-ray emission, strong variability, or X-ray over- 
luminosity) that could hint at binarity or the effects of magnetic 
confinement of the stellar wind. 



4.3. MT91516 

This star is of spectral type 05.5 V((f)). Its fundamental param- 
eters were determined by Herrero et al. (1999). Scaled to the 
distance adopted in the present paper, the bolometric luminos- 
ity is equal to log ^ = 6.01. Proceeding in the same way with 
the bolometric magnitude inferred by Negueruela et al. (2008 ) 
yields log ^ =6.12. 
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G. Rauw: A multi-epoch XMM-Newton campaign on Cyg OB2 
Table 8. Best fits of the EPIC spectra of Cyg OB2 #12 for photon energies above 0.4keV. 



CygOB2 #12, wabs x (apec(T,) + apec(T 2 )) 



Obs. 


Inst. 


d.o.f. 


x 1 


N H 








norm! 




kT 2 




norm 2 
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Jx 
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Jx 










(10 22 cm- 2 ) 


(keV) 


(10- 14 cm" 5 ) 




(keV) 


(10-' 4 cm- 5 ) 




(erg 


zm 1 s~' 




(ergcirr 2 s -1 ) 


1 


Ml+2 
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1.29 
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§) x 10- 
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3 
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) x 10- 
f) x 10- 


12 


2.64 x 10-" 


4 


EPIC 
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1.26 


1.83+ 


111 


0.74+ 


Hi 
Yo 


(8.84+° 


*}) x 10- 


-3 


1.81+ 


v;: 


3 


(2.74+ 


12 


2.61 x 10-" 


5 


EPIC 


108 


1.15 


2.06+ 


lis 


0.78 + 


(7.32+1 


|) x 10" 


-3 


2.28+ 


n 


3 


(1.99+ 


f ) x 10- 


12 


2.00 x 10-" 


6 


EPIC 
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II? 

08 


0.72 + 


IS 
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-3 
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31 
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12 
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Notes. The first three columns yield the number of the pointing, the instrument combination used in the fit and the number of degrees of freedom. 
The reduced,^ 2 is given in column 4, whilst the columns 5 to 9 provide the parameters of the absorption component and the emitting plasma. The 
last two columns yield the observed and absorption corrected flux over the 0.5 - lO.OkeV energy range. 




Time (JD-2450000) 

Fig. 9. Left: EPIC spectra (along with the best-fit model) of CygOB2 #12 as observed during the third observation. Right: Best-fit 
parameters of the EPIC spectra of Cyg OB2 #12. The top panel illustrates the observed flux (in units lCT 12 erg cirr 2 s -1 ) over the 
0.5 - 10 keV energy range. The middle panel yields the column density towards the X-ray emitting plasma, whilst the bottom panel 
displays the two temperatures of the emission components. 



From the colour excess E(B - V) = 2.48 (Massey & 
Thompson |1991| l, we infer a large interstellar column density of 
A^ SM = 1.44 x 10 22 cm- 2 . Nevertheless, the spectra display evi- 
dence of additional absorption that we modelled by including an 
ionized wind component. To achieve a reasonably good fit to the 
EPIC spectra, a two -temperature plasma model is required. The 
results of the fit to the EPIC spectra of MT91 516 (at energies 
above 0.4 keV) are listed in TablefTT] 

There is a degeneracy between the value of the ionized-wind 
column density and the temperature of the softer plasma com- 
ponent, and our best-fit solution oscillates between kT\ near 
0.7 keV with a very low wind column density (observations 2 
and 4), and kTi near 0.25 keV with a significantly larger wind 
column density (all other observations). The second temperature 
component at kT2 = (2.0 + 0.2) keV is more tightly constrained. 

The MOS 1 and pn count rates of this source are found to be 
variable though with a moderate amplitude. The same kind of 
variability is found in the observed flux of the source, which dis- 
plays variations with an amplitude of about 22% (peak to peak). 



Owing to the degeneracy of kT\, the ISM absorption-corrected 
flux displays very large variations. The mean value of the ISM 
absorption-corrected X-ray flux (3.71 ± 1.84) 10 -12 ergcirr 2 s -1 
yields log ^ = -6.79 ± 0.22 (for the Herrero et al.QjSDbolo- 
metric luminosity). Within the (large) errorbars, we find no sig- 
nificant X-ray overluminosity, relative to the canonical j^- value 
(Naze 120091 . 

The star displays radial velocity variations with a peak-to- 
peak amplitude of 113kms-' (Kiminki et al. [2007 ). Although 
these variations have an 11.5% probability of being caused by 
random noise, they might also represent binarity. Mason et al. 
(2009) reported the detection of an astrometric companion of 
similar brightness (magnitude difference of 0.4) at a separation 
of 0.7 arcsec. There is currently no evidence of any physical as- 
sociation between the two stars. Since there was no change in the 
position angle and separation of this pair over the time-frame of 
our X-ray observing campaign (Mason et al. 2009), it is unlikely 
that the observed X-ray variability, and by virtue the X-ray emis- 
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Table 10. Best-fits of the EPIC spectra of Cyg OB2 #7 for photon energies above 0.4 keV. 



Obs. Inst. d.o.f. 'x 1 



Cyg OB2 #7, wabs x wind x apec(T) 
logif 53 



fobs 

Jx 

(ergcrrT 2 s -1 ) 
(1.73+ | 3 ) x fir 13 
(1.78! 2 ° )x 10- 13 
(1.70+*) x 10- 13 
(1.75!;'°)x lO" 13 
(1.96^ 8 )xl0- 13 
(1.66 + ^ 3 ) x lO" 13 



JX 

(ergcnT 2 s~') 
8.81 x 10" 13 
7.53 x lO -13 
9.17 x 10 13 
10.97 x 10" 13 
7.60 x 10- 13 
8.49 x 10~ 13 



1 


EPIC 


59 


1.03 


21 


2 


Ml+2 


34 


1.25 


21 


3 


EPIC 


73 


1.75 


21 


4 


Ml+2 


22 


0.71 


21 


5 


Ml+2 


20 


0.88 


21 


6 


Ml+2 


28 


1.14 


21 



( cm' 2 ) 

1« 

• /u -.08 
7Q+.05 
• ' y -06 

7 2+.i7 

.83 + jo 
84 +:m 
,<¥+ -.io 



kT 
(keV) 

76 + 07 
0.60!* 
0.56! : °» 
77 +:i6 

'-.10 

0.58!;JI 



norm 
(1Q- I4 cnr 5 ) 

(1.27^) xlO- 3 
(0.89!° -\l)x 1Q- 3 
(1.30!° : §) x 10- 3 
(1.39!^) x lO" 3 
(1.00!° 32 )x 10- 3 
(1.40+^) x 10" 3 



Notes. The fit was perfomed with N^ M fixed to 1 .Ox 1 22 cm 2 . The first three columns yield the number of the pointing, the instrument combination 
used in the fit, and the number of degrees of freedom. The reduced x 1 is given in column 4, whilst the columns 5 to 9 provide the parameters of 
the absorption component and the emitting plasma. The last two columns yield the observed and ISM absorption corrected fluxes over the 0.5 - 
10.0 keV energy range. 



Table 11. Same as Table 10 but for MT91 516. 



log fr™ kT 1 norm! kT 2 norm 2 f° bs fi 



MT91 516, wabs x wind x (apec^O + apec(T 2 )) 



Obs. 


Inst. 


d.o.f. 


x 1 


1 


EPIC 


57 


1.18 


2 


EPIC 


82 


0.84 


3 


EPIC 


52 


1.2 


4 


EPIC 


53 


1.09 


5 


EPIC 


20 


1.13 


6 


EPIC 


35 


0.88 



X 



(cm- 2 ) (keV) (lQ-' 4 cnr 5 ) (keV) (lQ-' 4 crrr 5 ) (erg cirr 2 s' 1 ) (ergcnr 2 s-') 

21.5! 3 0.28!^ (5.72^ s 72 4 ) x 10~ 3 2.22! 3 1 ' (6.32^-$) x 10~ 4 (4.22!^) x 10~ 13 4.55 x 10~ 12 

19.6 +1 ^ 0.77 +:22 (6.98 +2 ^ 3 ) x lO" 4 2.26!^ (6.00 + ° : ^ 4 ) x lO" 4 (4.97!° ° 3 ) x lO" 13 2.48 x lO" 12 

21.7!' 2 0.22+ : '° (14.5™) x 10~ 3 1-91! 32 (8.85!™) x 10" 4 (4.58!^) x lO" 13 7.04 xlO" 12 



Notes. The fit was perfomed with N^ M fixed to 1.44 x 10 22 cm" 



19.0! 2 ^ 0.61! : >° (8.19!° p)x 10~ 4 1.86! 3 ° (7.25!j : = 3 ) x 10~ 4 (4.73!™) x 10" 13 2.97 x lO" 12 
22.1! 3 0.27! :3 = (18.0!ff 4 °)x lO" 3 1.85!™ (8.99!™) x 10~ 4 (4.14!™) x 10~ 13 2.08 x lO" 12 
21.7 +:2 0.32! : ;| (4.84 +33 4 « 1 )x 10- 3 1.90+^ (6.72 +2 |«) x 10~ 4 (3.96 + °° 3 ) x lO" 13 3.14 xlO" 12 



sion itself, could be associated with interactions between the as- 
tro me trie components. 



4.4. The pair Cyg OB2 #22 A + B 

CygOB2 #22 was classified as an 04111(f) giant and a fit to its 
spectrum yielded log — 6.04 (from the results of Herrero 
et al. 119991 scaled to our distance). The bolometric magni- 
tude given by Negueruela et al. (2008 ) yields a slightly larger 
log ^ = 6.10. However, CygOB2 #22 was subsequently re- 
ported to be a visual binary system consisting of an 03 1 (star A) 
and an 06 V star (star B, Walborn et al. l2002l Mason et al. l2009l 
separated by about 1.5arcsec. The 06 V component is itself a 
double system with a separation of 0.2 arcsec and a magnitude 
difference of 2.34 in the z filter (Sota et al. 1201 11 >. A priori, this 
result casts some doubt on the luminosity inferred by Herrero et 
al. ([1999 ), although the sum of the 'typical' luminosities of an 
03 1 and an 06 V star (Martins et al. 2002 ) amounts to a very 
similar value log ^ = 6.07. In the following, we adopt this 
latter value. We note that the separation between the A and B 
components is too small to resolve their X-ray emission with 
XMM-Newton. 

From the B - V — 2.04 colour index (Massey & Thompson 
119911 1. we derive an interstellar column of A^ SM = 1.35 X 
10 22 cm 2 . To achieve a good fit of the EPIC spectra, we need to 
include an ionized wind absorption. In addition, a second plasma 
temperature is required for the first four observations. In the last 
two observations, which have lower fluxes and lower signal-to- 



noise ratio, this second temperature is not constrained and is ac- 
tually not required to achieve an acceptable fit. 

The observed flux varies by nearly a factor of two between 
the second and the sixth observations. This variability is not due 
to an instrumental effect as it is also seen in the count rates of all 
the individual EPIC cameras. The EPIC count rates of this source 
are indeed found to be variable at the 99.9% level (though not 
in all energy bands for each instrument). The ISM absorption- 
corrected fluxes vary by a factor of three and the extreme values 
correspond to Lx in the range from 1 .5 x 10 32 to 4.5 x 10 32 erg s _1 
and log 7^ between -7.47 and -7.00. The star hence varies from 
underluminous by a factor three to just at the typical X-ray lu- 
minosity level. 

4.5. CPR2002A11 

Negueruela et al. (2008 ) inferred log ^ = 5.44 (again scaled 
to the distance adopted in the present paper) for this 07.51b- 
11(f) star. From B - V = 2.19 (for MT91 267), we estimate that 
N^ M = 1.43 x 10 22 cirr 2 . The spectra are unlikely to contain 
an additional (wind) absorption component, thus we fitted them 
using a single-temperature plasma model with fixed interstellar 
absorption. 

This source is found to be highly variable in count rate 
(in all three EPIC instruments) and this variability is reflected 
in the observed flux that varies by a factor of 2.7. The ISM 
absorption-corrected fluxes yield log j^- between -6.84 and 
-6.44. This star hence appears somewhat overluminous in X- 
rays. Furthermore, that its X-ray emission varies and is domi- 
nated by a relatively hot plasma at (1.60 ± 0.15)keV suggests 
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Table 12. Same as TablefTol but for Cyg OB2 #22. 



logA ,wmd 

(cm" 2 ) 



CygOB2#22, wabs x wind x (apec(T,) + apec(T 2 )) 



Obs. 



Inst. 



d.o.f. 



1 


EPIC 


62 


1.18 


2 


EPIC 


56 


0.95 


3 


EPIC 


70 


1.18 


4 


EPIC 


37 


1.69 


5 


M1+M2 


9 


1.45 


6 


EPIC 


22 


1.23 



21.6+' 
21.6 +0 3 1 

21.8!° 3 2 
21.7+ 3 
21.53 



kT, 
(keV) 
0.50+ 



norm[ 
(lQ-' 4 cm- 5 ) 
) x IO" 3 
')x 10- 



5+3155T 



kT 2 
(keV) 



2.27 
3.33 
0.89 
1.62 



■88 

.10 

-.i i 

-.85 
-.38 



norm 2 

(10-' 4 cm- 5 ) 



fobs 
JX 



(2.10 +1 8 28 )x 10- 4 
(2.52^«) x 10- 4 
(6.96! 2 g) x 10- 4 
(4.00+ 2 * 5 ) x 10- 4 



(erg cm 2 s 1 



(2.64!" 2,) ) x HF 3 
(3.48S)xl0- 13 
(l-99^5 2 )x 10 13 
(2.41^ : ;0)x io- 13 
(1.96*™) xlO- 13 
(1.77^) x 10-' 3 



run 
JX 

(erg cirr 2 s -1 



1.99 x IO" 12 
1.94 x IO" 12 
2.15 x IO-' 2 
2.62 x IO" 12 
0.88 x IO" 12 
1.31 x IO" 12 



0.57 
0.30; 



:/>?> 

.10 
.10 

0.29+° 4 
0.93! :32 
0.73! :32 



(1 
(1 

( 3 - 16 • 



(1.5l ;t) - 7( h- ,n-3 



') X 10- 



(6.51^) x 10- 3 
(0.91 + ^ 



;■«) x 10- 

■'M\ ~ 1 n- 



(1.05!*j)x 10 



Notes. The fit was perfomed with N™ M fixed to 1.35 x 10 22 cm- 



Table 13. Same as Table 10 but for CPR2002 Al 1 



CPR2002A11, wabs x apec(T) 



Obs. 


Inst. 


d.o.f. 


x 1 


kT 


norm 




fobs 
JX 


run 
JX 












(keV) 


(io- 


14 cm- 5 ) 




(ergctrr 2 s -1 ) 


(erg cm -2 s 


"') 


1 


EPIC 


85 


1.60 


1 63 +u? 

-06 


(1.20+ 




3 


(5.29+ : {\) x IO" 13 


1.70 x 10" 


12 


2 


EPIC 


52 


1.27 


I 49+.M 


(0.60+ 




-3 


(2.44 +: !«) x IO" 13 


0.88 x 10" 


12 


3 


M1+M2 


51 


1.89 


1 52 +: ^ 
-•98 


(1.19+ 


5p x 10 


-3 


(4.96 + 2 «) x IO" 13 


1.74 x 10" 


12 


4 


M1+M2 


25 


1.67 


1.89^ 


(1.05 + 


)X10- 


-3 


(5.21! :43 )x 10- 13 


1.43 x 10" 


12 


5 


M1+M2 


24 


1.07 


1 57 + - IS 
-44 

i.49!;= 


(1.56 + 


») x 10 
*f)xlO- 


-3 


(6.62+ |h x IO" 13 


2.22 x 10" 


12 


6 


M1+M2 


22 


0.94 


(0.97+ 


-3 


(3.95 + 33 ) x IO" 13 


1.42 x 10" 


12 



Notes. The fit was perfomed with A^ SM fixed to 1.43 x 10 22 cm" 2 . 



that it might either be an interacting wind binary system or a star 
containing a magnetically confined wind. Spectroscopic moni- 
toring in the optical is highly desirable to help us establish the 
multiplicity of this object. 



5. Conclusions 

We have discussed the X-ray properties of a number of stars in 
the core of Cyg OB2. One of this association's most interesting 
stars is the LBV candidate Cyg OB2 #12. We have shown for the 
first time that the X-ray flux of this star varies on many differ- 
ent timescales ranging from days to years. These variations are 
mainly due to variations in the column density and changes in 
the emission measure of the harder component of the spectrum 
as one would expect for models of either an eccentric colliding 
wind binary or a rotationally modulated magnetically confined 
wind. Given that the variability occurs on timescales as short as 
days, the latter scenario seems slightly more plausible. An ex- 
tensive monitoring of this star, both in the optical and the X-ray 
domain could shed new light on the nature of this intriguing ob- 
ject. 

Apart from the known colliding-wind systems, the majority 
of the O-type stars in CygOB2 appear to have rather normal 
X-ray properties. Their comparatively hard X-ray emission is 
not intrinsic to the sources, but is caused instead by the huge 
interstellar absorption towards the association. Several objects 
feature a second thermal component that is too hot to be easily 
explained by the wind-embedded shocks scenario. However, it 
must be stressed that this result is not restricted to the stars in our 
sample, nor direct evidence of either multiplicity or a magnetic 
field. Indeed, a second thermal component with a temperature 
around 2 ke V is frequently required to fit high quality CCD X- 
ray spectra of O-type stars, regardless of their multiplicity (e.g. 



Naze et al. 1201 11 1. In our sample, the exceptions to this general 
statement are MT91 516 and CPR2002 All, which we propose 
to be good candidates for binary systems. About one third of the 
OB stars display variability between our pointings. All known 
binary systems (except for the W UMa system Cyg OB2 #27) are 
highly variable, which suggests that the other variable OB stars 
might also be binaries. However, we stress that these objects do 
not necessarily display a permanent excess X-ray emission; the 
Trapezium-like system Cyg OB2 #22 even appears to be 'X-ray 
underluminous' in some of our observations. 

Finally, a large number of the secondary sources, the bulk of 
which are most likely associated with pre-main sequence low- 
mass stars, also display highly significant variability. However, 
this variability appears to be more complex than the impulsive 
rise followed by an exponential decay that is commonly seen in 
flaring stars. Clarifying the relationship between these objects 
and the massive star population in Cyg OB2 will help us deter- 
mine the impact of massive stars on star formation in their sur- 
roundings. In this respect, the Chandra very large program sur- 
vey of Cyg OB2 (Drake [201 \\ and its multi-wavelength support 
observations will certainly help us make substantial progress. 
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